function result=Fig1_f_1()
filename = 'gbsweep_152.GBs';
pathname = 'raw data/';

[data, ~] = load_GBs(filename, pathname);

data.gate = data.gate - (-20e-3); %new Vb = original Vb - Vm
n_bias = length(data.bias);
n_gate = length(data.gate);
peak = zeros(1,n_gate);
width = zeros(1,n_gate);
result.gate_2 = [];%pure
result.peak_2 = [];
result.gate_1 = [];%remove bad fit
result.peak_1 = [];
result.valid = [];
width_2 = [];
result.nu = [];

for i=1:n_gate
    try
        [~,tmp_max] = max(data.map(:,i));
        tmp = fit(data.bias(tmp_max-3:tmp_max+3)',data.map(tmp_max-3:tmp_max+3,i),'gauss1');
    catch
        tmp.b1 = 0; tmp.c1 = 1;
    end
    width(i) = tmp.c1;
    peak(i) = tmp.b1;
    
    result.gate_2 = [result.gate_2 data.gate(i)];
    result.peak_2 = [result.peak_2 peak(i)];
    
    %remove peaks that is wider than 10mV (the peak is out of range)
    if width(i)< 10e-3 
        result.valid = [result.valid 1];
        result.gate_1 = [result.gate_1 data.gate(i)];
        result.nu = [result.nu -0.5 + (data.gate(i)+1.4025)/0.555 ];
        result.peak_1 = [result.peak_1 peak(i)];
    else
        result.valid = [result.valid 0];
    end
    
    
end

Fcolor='k';
Bcolor='w';
Fsize=18;
Bsize=14;

f=figure();
plot(result.nu,-1*(1e3*(result.peak_1)-27.37),'LineWidth',2);

f.Color=Bcolor;
ax = gca;
ax.FontSize=Fsize;
ax.Color = 'none';
ax.XColor=Fcolor;
ax.YColor=Fcolor;
yl=ylabel('u (meV)');
xl=xlabel('nu');
yl.FontSize=Fsize;
xl.FontSize=Fsize;

x_limit=xlim;
x_limit=[-0.85,-0.15];
y_limit=ylim;
y_limit=[-5,5];
axis([x_limit(1),x_limit(2),y_limit(1),y_limit(2)]);

end